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Abstract 

A  multi-physics  model  coupling  electrochemical  kinetics  with  fluid  dynamics  has  been  developed  to  simulate  the  transport  phenomena  in 
mono-block-layer  built  (MOLB)  solid  oxide  fuel  cells  (SOFC).  A  typical  MOLB  module  is  composed  of  trapezoidal  flow  channels,  corrugated 
positive  electrode-electrolyte-negative  electrode  (PEN)  plates,  and  planar  inter-connecters.  The  control  volume-based  finite  difference  method 
is  employed  for  calculation,  which  is  based  on  the  conservation  of  mass,  momentum,  energy,  species,  and  electric  charge.  In  the  porous 
electrodes,  the  flow  momentum  is  governed  by  a  Darcy  model  with  constant  porosity  and  permeability.  The  diffusion  of  reactants  follows 
the  Bruggman  model.  The  chemistry  within  the  plates  is  described  via  surface  reactions  with  a  fixed  surface-to-volume  ratio,  tortuosity  and 
average  pore  size.  Species  transports  as  well  as  the  local  variations  of  electrochemical  characteristics,  such  as  overpotential  and  current  density 
distributions  in  the  electrodes  of  an  MOLB  SOFC,  are  discussed  in  detail. 

©  2004  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Fuel  cells  are  drawing  keen  attention  for  their  high  ef¬ 
ficiency  and  eco-friendly  characteristics.  Among  fuel  cells, 
solid  oxide  fuel  cells  (SOFCs)  are  less  sensitive  to  fuel  com¬ 
position  compared  to  other  fuel  cell  systems.  Besides  the 
reformed  gas  from  the  external  reformer,  the  anode  feed  of 
SOFC  can  employ  natural  gas,  LPG,  methanol,  or  coal  gas 
through  internal  reforming  on  the  anode  side. 

A  typical  SOFC  operates  at  a  temperature  between 
700  and  1000  °C.  Such  a  high  temperature  accelerates  the 
electrochemical  reaction;  therefore,  SOFCs  do  not  require 
precious  metal  catalysts  to  enhance  the  reaction.  In  addi¬ 
tion,  SOFCs  can  attain  an  electric  efficiency  as  high  as 
65%,  surpassing  the  50%  achieved  by  LNG  compound- 
power  generation  because  of  its  ability  to  be  used  in  com¬ 
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bination  with  gas  turbines  using  high-temperature  gas  ex¬ 
haust. 

SOFCs,  like  any  other  electrochemical  device,  consist 
of  two  electrodes  (anode  and  cathode),  and  an  electrolyte, 
which  allows  transfer  of  oxygen  ions.  These  three  layers  are 
laminated  into  a  single  positive  electrode/electrolyte/negative 
electrode  (PEN)  plate.  The  PEN  plates  are  then  connected  in 
a  series  stack  via  inter-connect  plates.  As  shown  in  Fig.  1,  the 
PEN  plate  can  be  molded  unevenly  in  a  trapezoidal  or  dimple 
shape,  which  is  known  as  mono-block  layer  built  (MOLB) 
SOFC.  Such  a  design  has  benefits  in  enhancing  its  mechani¬ 
cal  strength  and  increasing  its  power  density.  The  corrugated¬ 
shaped  PEN  plate  not  only  ensures  that  the  effective  active 
area  is  higher  than  the  projected  area,  but  also  provides  the 
film  with  the  combined  function  of  fuel  and  airflow  paths, 
making  the  cell  compact. 

For  optimal  design  and  operation  of  MOLB-type  SOFCs,  a 
comprehensive  understanding  of  transport  and  electrochem¬ 
ical  characteristics  is  highly  essential.  However,  previous 
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Nomenclature 

a'  stoichiometric  coefficient  of  the  products 

a"  stoichiometric  coefficient  of  the  reactants 

D  diffusivity  (m2  s-1) 

F  Faraday’s  constant  (96487  C  mol-1) 

i  current  density  (A  m-2) 

j  transfer  current  density  (A  m-3) 

J  diffusive  flux  (kg  s-1  m-2) 

k  thermal  conductivity  (W  K- 1  m- 1 ) 

M  molecular  weight  (kg  mol - 1 ) 
p  pressure  (Pa) 

R  universal  gas  constant  (W  mol-1  K-1) 

T  temperature  (K) 

U  velocity  vectors  (ms-1) 

X,  y,  Z  coordinate  system,  Fig.  1  (m) 

Xi  fraction  of  molar  concentration  of  species  i 
Yi  mass  fraction  of  species  of  i 

Greek  symbols 

0  solid-phase  potential  (V) 

a  symmetric  factor 

o  electric  conductivity  ( Q ~ 1  m- 1 ) 

s  porosity 

k  permeability  (m2) 

p  density  (kg  m- 3 ) 

0  (or  Over_P)  overpotential,  0q  —  (V) 

li  viscosity  (m  s-2) 

Subscripts 
A  anode 

C  catalyst  or  cathode 

eff  effective 

E  electrolyte 

F  fluid 

i  species 

ref  reference 

S  solid 

T  transfer  current 


studies  are  mostly  for  either  planar  or  tubular  SOFCs  [1-7], 
and  so  far  have  not  yet  considered  the  geometry  of  MOLB- 
type  SOFCs.  The  present  study  aims  to  develop  a  CFD-based 
model  to  resolve  the  complex  interaction  between  fluid  dy¬ 
namics  and  electrochemical  kinetics  in  a  MOLB-type  SOFC. 
In  contrast  to  the  results  obtained  from  models  that  specify 
uniform  surface  overpotential  [4,8],  the  present  model  im¬ 
plements  the  voltage-to-current  algorithm  that  allows  for  a 
more  realistic  spatial  variation  of  electrochemical  kinetics. 
This  results  in  an  improved  prediction  of  local  overpotential 
and  current  density  distributions. 


2.  Model  development 

2.7.  Basics  of  MOLB-type  SOFC 

A  typical  single-cell  module  of  the  MOLB-type  SOFC 
consists  of  a  PEN  plate  and  an  inter-connect  (Fig.  1).  The 
fuel  and  air  run  in  a  counter  current  direction  in  the  trape¬ 
zoidal  channels  that  are  separated  by  a  corrugate  PEN  plate. 
The  oxidant  reduction  reaction  occurring  in  the  cathode  is 
expressed  as  follows. 

O2  +  4e“  — 202-  (1) 

The  oxygen  ion  is  then  transported  through  the  electrolyte 
and  then  into  the  active  anode.  The  electrochemical  reactions 
of  H2  and  CO  on  anode  catalysts  are: 

H2  +  02“  -*  H20  +  2e“  (2) 

CO  +  O2-  -*  C02  +  2e“  (3) 

The  co-oxidations  of  Eqs.  (2)  and  (3)  are  limited  to  a 
single-step  reaction  in  the  present  model,  i.e., 

H2  +  CO  +  202-  -*  H20  +  C02  +  4e”  (4) 


2.2.  Model  equations 


The  continuity,  momentum,  and  species  conservations  are 
modeled  as: 


V  •  (epU)  =  0  (5) 

e2nU 

V  •  (epUU)  =  —sVp  +  V  •  (sfieffVU)  H - — —  (6) 

K 


V  •  (epUYi)  =  V  •  J i  +  (a"  -  a'jf  (7) 

The  isotropic  linear  resistance  model  of  Darcy’s  law  with 
constant  porosity  (s)  and  permeability  ( k )  is  used  to  govern 
the  momentum  equation.  The  diffusive  mass  flux  J/  is  de¬ 
scribed  by  the  Stefan-Maxwell  equation  [9].  The  effective 
diffusivities  in  the  multi-component  mixture  system  follow 
the  Bruggman  model  [10].  a"i  and  a'i  are  the  stoichiometric 
coefficients  of  the  products  and  reactants,  respectively.  The 
electrode  electrochemistry  and  the  species  transport  kinet¬ 
ics  are  connected  by  the  Butler- Volmer  equation  [1 1]  on  the 
surfaces  of  the  active  layers,  i.e., 


7t,a  —  7o,a  ' 
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Fig.  1.  Schematic  diagram  of  a  single-cell  module  of  the  MOLB-type  SOFC. 
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where  r\  is  the  electrode  overpotential,  jo,a  and  yo,c  the  an¬ 
odic  and  cathodic  exchange  current  densities,  respectively, 
and  a  a  and  ac  are  the  anodic  and  cathodic  Tafel  constants, 
respectively. 

In  the  active  layer,  a  potential  difference  exists  between 
the  catalyst  (<Pq)  and  electrolyte  (0E)  to  drive  the  transfer 
current  (/t,a  or  jt,c)>  keeping  the  electrochemical  reaction 
continuously.  As  shown  in  Fig.  2,  the  current  passes  through 
the  electrode  can  be  decomposed  two  parts,  i.e., 


i  =  ic  +  *e 


(10) 


These  two  current  components  interact  through  the  elec¬ 
trochemical  reaction.  The  electrons  are  either  transferred  to 
the  catalyst  from  the  electrolyte  or  vice  versa  during  the  reac¬ 
tion.  Application  of  Ohm’s  law  to  Eq.  (10)  yields  the  current 
conservation: 


V  •  (ic)  =  V(-aeff5cV0c)  =  ~j r 


(11) 


V  •  (*E)  =  V(-aeff,EV0E)  =  jT  (12) 

where  <reff,c  and  aeff,E  are  the  effective  electric  conductivities 
of  the  catalyst  and  electrolyte  in  the  electrode,  respectively. 
Physical  properties  and  parameters  used  in  the  present  model 
are  summarized  in  Tables  1  and  2. 


3.  Numerical  implementation 

As  shown  in  Fig.  3,  a  typical  MOLB  unit  is  symmetric 
about  the  mid-planes  of  the  air  and  fuel  channels.  There¬ 
fore,  the  present  computational  domain  covers  the  region  be¬ 
tween  these  two  planes  for  memory  and  time  saving.  The 
governing  equations  were  numerically  solved  by  the  control 
volume-based  finite  difference  method  [12].  The  discretiza¬ 
tion  procedure  ensures  conservation  of  mass,  momentum,  and 
concentration  over  each  control  volume.  The  coupling  of  the 
pressure  and  velocity  fields  was  treated  via  the  SIMPLER 
pressure  correction  algorithm  [13,14].  The  upwind  differ¬ 
ence  scheme  was  used  to  treat  the  convection  and  diffusion 


Ionic  Transfer 


density,  /E  density,  /c 


Fig.  2.  Scheme  of  the  transfer  current  in  the  SOFC  electrode. 
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Table  1 

Porous-chemistry  properties  of  the  electrodes 


Description  Unit  Value 


Anode  exchange  current,  /o,a 

Am  3 

1.0  X  109 

Cathode  exchange  current,  /o,c 

Am-3 

5.0  x  107 

Anode  transfer  coefficient,  a  a 

— 

1.0 

Cathode  transfer  coefficient,  ac 

— 

1.0 

Anode  stoichiometric  coefficient 

— 

3.0 

Cathode  stoichiometric  coefficient 

— 

5.0 

Anode  inlet  mass  flow  rate 

kgs-1 

1.19  x  10“6 

Cathode  inlet  pressure 

kgs-1 

4.74  x  10“8 

Anode  inlet  pressure 

kPa 

1.013  x  105 

Cathode  inlet  pressure 

kPa 

1.013  x  105 

Stack  temperature 

K 

973 

Species  molar  concentration  at  anode  inlet 

Hydrogen,  Xn, 

% 

45.57 

Carbon  monoxide,  Xco 

% 

6.97 

Carbon  dioxide,  Xqo2 

% 

6.55 

Water,  Xh2o 

% 

36.84 

Methane,  Xch4 

% 

4.07 

Total 

% 

100 

Species  molar  concentration  at  cathode  inlet 

Oxygen,  X02 

% 

21.11 

Nitrogen,  X^2 

% 

78.99 

Total 

% 

100 

terms.  Solutions  were  considered  to  be  converged  after  the 
ratio  of  residual  source  (including  mass,  momentum,  species, 
and  current)  to  the  maximum  flux  across  a  control  surface  is 
less  than  1.0  x  10-3.  The  computations  were  performed  on 
18  x  40  x  50  (X  by  Y  by  Z)  structured,  orthogonal  meshes 
in  the  present  work  (Fig.  3).  A  typical  simulation  required 
about  500  min  of  central  processing  unit  time  on  a  Pentium 
IV  2.0  GHz  PC. 

4.  Results  and  discussion 

Conditions  for  the  present  calculation  are  to  achieve 
average  cell  temperature  of  973  °C.  Air  and  fuel  deliv¬ 
ered  to  the  cathode  and  anode  are  2.975  x  10  5  kgs  1 
and  1.185  x  10  6 kgs  1 ,  respectively.  These  correspond 
to  the  stoichiometric  flow  ratios  5.0  and  3.0  at  the  current 
density  of  1.0  A  cm-2,  respectively.  The  molar  composi¬ 
tion  of  the  fuel  is  45.57%/6.97%/6.55%/36.84%/4.07%  for 
H2/CO/CO2/H2O/CH4.  The  voltage  drop  across  the  module 
is  set  at  0.4  V  in  the  present  computation. 


Fig.  3.  Dimensions,  coordinate  system,  and  grid  distribution  of  the  compu¬ 
tational  domain. 


Fig.  4  displays  the  distribution  of  oxygen  mass  fraction 
(To2)  on  the  interface  between  the  air  channel  and  the  cath¬ 
ode.  It  is  observed  that  the  high  oxygen  mass  fraction  from 
the  entrance  of  the  channel  inlet  decreases  along  the  axial 
distance  (— Z)  due  to  the  oxygen  reduction  reaction.  The  oxy¬ 
gen  mass  fraction  has  a  local  minimum  at  the  channel  corner 
formed  by  the  inclined  and  bottom  planes.  The  forced  con¬ 
vection  in  the  local  minimum  is  weak  due  to  the  corner  effect. 
Another  corner  formed  by  the  inclined  plane  and  the  upper 
plane  does  not  sense  the  local  minimum.  It  is  clearly  seen 


Table  2 


Physical  properties  in  the  porous  electrodes 


Description 

Anode 

Cathode 

Electrolyte 

Electrode  contact 

Inter-connector 

Thickness  (mm) 

1.4 

0.4 

0.2 

0.4 

0.5 

Porosity  s 

0.4 

0.5 

1.0  X  10“5 

0.3 

1.0  X  10“5 

Permeability,  k  (m2) 

1.0  X  10“12 

1.0  X  10“12 

1.0  x  10“18 

1.0  X  10"12 

1.0  x  10“18 

Tortuosity 

1.5 

1.5 

6 

1.5 

6 

S/VCm”1) 

1000 

1000 

— 

— 

— 

Pore  size  (m) 

1.0  X  10“6 

1.0  X  10“6 

1.0  x  10“6 

1.0  X  10"6 

1.0  x  10“6 

Effective  thermal  conductivity,  keff  (W  m~ 1 ) 

6.23 

9.6 

2.7 

60.7 

9.6 

Effective  electric  conductivity,  (7eff  (£2 1  m-1) 

1.0  X  105 

7.7  X  103 

1.0  X  10“2° 

1.56  x  107 

1.35  X  106 
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Fig.  4.  Oxygen  concentration  distributions  at  the  interface  of  the  channel  and  the  cathode. 


Fig.  5.  Reactant  concentration  distributions  on  the  interface  between  channel  and  the  anode  (a)  hydrogen  and  (b)  carbon  monoxide. 
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that  the  iso-line  of  O2  mass  fraction  extends  from  the  mid- 
inclined-plane  to  the  plateau  center,  which  is  in  touch  with  the 
cathode  contact  through  the  upper  channel  corner.  Since  the 
air-flow  is  retarded  by  the  porous  matrix  of  the  cathode  con¬ 
tact,  the  fresh  air  cannot  move  into  the  upper  plateau  easily. 
This  results  in  a  lack  of  reactants  for  electrochemical  reaction. 
In  this  circumstance,  a  significant  overpotential  will  arise  to 
operate  normally  in  this  region,  which  will  be  shown  later. 
It  can  be  concluded  that  the  oxygen  transport  is  dominated 
by  convection  on  both  surfaces  of  the  inclined  and  bottom 
planes.  In  contrast,  the  oxygen  transport  along  the  plateau  is 
mainly  achieved  by  diffusion. 

Fig.  5  shows  the  distribution  of  mass  fraction  of  hydro¬ 
gen  (Fh7)  and  carbon  monoxide  (Fco)  on  the  surfaces  of  the 
channel-anode  interface.  It  can  be  seen  that  both  Fh2  and 
Fco  gradually  decrease  due  to  the  chemical  reaction  along 
the  channel  distance.  The  local  minima  of  Fh?  and  Fco  also 
occur  at  the  corner  formed  by  the  inclined  and  ceiling  planes 
and  at  the  mid-plane  of  the  bottom  plane  where  H2  and  CO 
are  transported  mainly  through  diffusion. 


Fig.  6a  and  b  shows  the  overpotential  distributions  on  the 
mid-planes  of  the  cathode  and  anode,  respectively.  In  the 
anode,  the  positive  overpotential  directs  the  transfer  current 
from  the  catalyst  to  the  electrolyte,  and  vice  versa  in  the  cath¬ 
ode.  In  general,  the  cathode  needs  higher  energy  in  order  to 
activate  the  reactants.  The  absolute  value  of  the  overpotential 
in  the  cathode  is  significantly  higher  than  that  in  the  anode. 
It  can  be  seen  that  the  overpotential  is  higher  near  the  exit 
of  the  oxidant  flow  where  it  coincides  with  the  regions  of 
lower  oxygen  concentration  on  the  cathode  plane.  Depletion 
of  reactants  results  in  an  increase  of  overpotential  along  the 
channel.  At  a  fixed  axial  distance,  the  maximum  overpoten¬ 
tial  occurs  at  the  plateau  center  where  the  oxygen  is  almost 
used  up  (Fig.  4). 

Fig.  7a  and  b  compares  the  current  density  distribution  on 
the  mid-planes  of  the  cathode  and  anode.  These  plots  cover 
only  the  mid-quarter  channel  results.  The  arrows  attached  on 
the  plane  surfaces  are  the  projections  of  the  current  density 
vectors  at  the  corresponding  position.  These  arrows  indicate 
not  only  the  direction  but  also  the  relative  magnitude  of  the 
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Fig.  6.  Overpotential  distributions  on  the  mid-planes  of  the  electrodes  (a)  cathode  and  (b)  anode. 
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Fig.  7.  Distributions  of  current  density  on  the  mid-plane  surfaces  of  electrodes  (a)  cathode  and  (b)  anode. 


currents.  Attention  is  placed  on  the  inclined  planes  first.  On 
the  inclined  plane  of  the  anode,  the  current  density  gradually 
decreases  along  the  current  direction,  since  the  current  trans¬ 
fers  from  the  high-potential  catalyst  (electronic  conductor)  to 
the  low-potential  electrolyte  (ionic  conductor).  Conversely, 
the  current  density  gradually  increases  along  the  current  di¬ 
rection  on  the  cathode  plane  because  of  the  current  accumu¬ 
lation  by  transferring  the  current  from  the  ionic  conductor 
to  the  electronic  conductor.  In  the  present  model,  the  anode- 
supported  PEN  plane  has  a  thick  anode,  thus,  resulting  in 
a  current  passage  with  a  larger  cross-sectional  area.  This  re¬ 
duces  the  current  density.  The  top  and  the  bottom  planes  have 
a  relatively  large  area  for  the  passing  through  current,  which 
directs  upward  for  both  electrodes.  As  a  result,  the  current 
density  is  rather  small  in  these  regions. 

5.  Conclusions  and  future  work 

This  paper  has  presented  a  computational  simulation  of 
the  fluid  dynamics  and  electrochemical  kinetics  in  a  MOLB- 


type  SOFC.  Significant  results  about  the  local  transport  char¬ 
acteristics  inside  the  MOLB-type  SOFC,  such  as  the  reactant 
gas  concentration  distribution,  the  overpotential  distribution, 
and  the  current  density  distribution  have  been  presented.  The 
unique  features  of  this  model  are  the  implementation  of  the 
voltage-to-current  algorithm  and  the  coupling  of  the  poten¬ 
tial  field  with  the  reactant  species  concentration  field,  which 
allows  for  a  more  realistic  spatial  variation  of  the  electro¬ 
chemical  kinetics.  The  spatial  variation  of  the  cathodic  and 
anodic  surface  overpotential  is  computed  locally,  thus,  re¬ 
sulting  in  an  improved  prediction  of  the  local  current  density 
distribution.  These  local  current  density  distribution  patterns 
are  radically  different  from  those  obtained  from  models  that 
do  not  account  for  the  non-uniformity  of  the  surface  overpo¬ 
tential. 
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